Crystal field, spin-orbit coupling and magnetism in a ferromagnet YTi03 
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Magnetic properties of stechiometric YTiOa has been calculated within the single-ion-based 
paradigm taking into account the low-symmetry crystal field and the intra-atomic spin-orbit cou- 
pling of the Ti 3+ ion. Despite of the very simplified approach the calculations reproduce perfectly 
the value of the magnetic moment and its direction as well as temperature dependence of the mag- 
netic susceptibility x(T). It turns out that the spin-orbit coupling is fundamentally important for 
3d magnetism and magnetic properties are determined by lattice distortions. 
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INTRODUCTION^ 

YTiC>3 is a unique 3d ferromagnet |]J, |2j, LI - the most 
of oxides are antiferromagnetic. In combination with 
a rather simple perovskite structure and Ti 3+ ions ex- 
pected to have one electron in the incomplete 3d shell, 
YTi03 is regarded to be very good examplary sys- 
tem for studying basic interactions in 3d oxides. De- 
spite this simplicity its properties are not understood 
yet. In fact, there is going on at present a hot debate 
on description of its magnetic and electronic properties 

A time when such system was treated as a S= 1/2 system 
i.e. with the spin-only magnetism and the fully quenched 
orbital magnetism is already gone. Also the simplest ver- 
sions of the band picture turned out to be completely use- 
less giving already at the start disagreement with exper- 
iment predicting a metallic ground state whereas YTiO,3 
is in the reality a good insulator. 

YT1O3, when stoichiometric, is very good insulator. 
Its resistivity at room temperature amounts to 5T0 -2 
ficm ^lj . The resistivity rapidly grows up with decreas- 
ingtemperature. It is ferromagnetic below T c of 30-35 K 
[20] . The macroscopic magnetisation, if recalculated per 
the formula unit, points to a moment of 0.84 ub 0- The 
paramagnetic susceptibility has been found to follow the 
Curie- Weiss law with 9 = 39 K, if substracted a diamag- 
netic and temperature independent orbital contributions 
of a total size of 0.351-10 -3 emu/mol. 

The aim of this paper is to present results of single-ion 
based calculations of properties of YTiC>3. Our single- 
ion results seem to be quite remarkable. We took into 
account low-symmetry off-octahedral crystal-field (CEF) 
interactions and the intra-atomic spin-orbit (s-o) cou- 
pling, that turn out to be of the comparable strength. 

THEORETICAL OUTLINE 

We consider exactly stechiometric YTi03. From this 
and the insulating ground state we infer that all Ti ions 
are in the trivalent state according to the charge distribu- 



tion Y 3+ Ti 3+ Og . The relevant charge transfer occurs 
during the formation of the compound. In the perovskite- 
based structure of YTi0 3 the Ti 3+ ion is surrounded by 
six oxygen ions forming distorted octahedron. There are 
still some Ti ions at the surface, for instance, with a 
reduced symmetry, but they are generally neglected, be- 
cause we are interested in intrinsic properties of YTiC>3. 
The local octahedra in YTiC>3 are tilted and rotated, 
what causes the need for consideration a larger elemen- 
tary cell, with four Ti ions instead of one as in the simple 
perovskite structure. Thanks it other three Ti ions get 
a crystallographic freedom. As a result of rotations, tilts 
and other atom displacements there are four short Ti-0 
bonds and two long Ti-0 bonds. 

The orthorhombic lattice in the Pbnm structure results 
from that of an ideal cubic perovskite by setting a =\/2 
a c , b a c and c =2a c , where a c denotes the lattice 

parameter of the simple cubic perovskite. a c is of order 
of 400 pm. 

In the Pnma structure used by some authors the dou- 
bling occurs along the b direction. The c, a, and b axis in 
the Pnma structure becomes the a, b and c axis in the 
Pbnm structure, respectively. Independently of the used 
Pbnm or Pnma space group the derivation of the local 
surroundings must be, of course, the same. 

Here we use the Pbnm space group. Then the easy 
magnetic axis is the c axis, The lattice parameters 
at T = 293 K, according to Ref. 2, cited by lEEl, are: 
a =531.6 pm, 6 =567.9 pm and c = 761.1 pm. These 
parameters have been confirmed by detailed structural 
measurements of Loa et al. 0] . 

Cwik et al. 0| have derived the respective bonds: 207.7 
pm (Ti-O(b)), 201.6 pm (Ti-0 fa)} and the apical bond 
202.3 pm (Ti-O(c)). Loa et al. |l9j have measured influ- 
ence of the external pressure on these bonds. From these 
crystallographic studies we get an input to our theoreti- 
cal considerations that the lattice surroundings of the Ti 
ion in YTi03 is predominantly octahedral with a slight 
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orthorhombic distortion. 



RESULTS AND DISCUSSION 

The electronic structure of the Ti 3+ ion with one d- 
electron under the action of the crystal-field interactions 
Hcf and in the presence of the intra-atomic spin-orbit 
coupling H s _ we calculated with the use of the single- 
ion-likc Hamiltonian often used in description of 3d im- 
purity states in the Electron Paramagnetic Resonance 
[22L l2.'l |24J, which we accept also for a solid, where 3d 
atom is the full part of the crystallographic structure: 

H d = H CF + H s - Q = 
Bi{Ol + 5(9|) + Xs-oL ■S + B%O 2 + Li B (L+g e S) ■ B ext .(l) 

The calculated electronic structure is presented in Fig. 
1. The crystal field has been divided into the cubic 
part, usually dominant in case of compounds containing 
3d ions, and the off-octahedral distortion written by the 
second-order leading term B^O". The last term, Zeeman 
term, allows calculations of the influence of the external 
magnetic field. g e amounts to 2.0023. The Zeeman term 
is necessary for calculations, for instance, of the param- 
agnetic susceptibility - in fact the paramagnetic suscepti- 
bility is customarily calculated [25] as the magnetization 
in an external field of, say, 0.1 T applied along different 
crystallographic directions. 

The detailed form of the Hamiltonian (1) is written 
down in the LS space that is the 10 dimensional spin- 
orbital space \LSL Z S Z ). The L and 5 quantum numbers 
for one 3d electron are equal to L—2 and 5=1/2 (here, for 
3d 1 configuration, lower I and s could be also used) . The 
Hamiltonian (1) is customarilly treated by perturbation 
methods owing to the weakness of the s-o coupling for the 
3d ions in comparison to the strength of the crystal-field 
interactions. We have accepted the weakness of the s-o 
coupling, what is reflected by the sequence of terms in 
the Hamiltonian (1), but we have performed direct cal- 
culations treating all shown terms in the Hamiltonian (1) 
on the same footing. The separate figures, if presented, 
are shown for the illustration reasons. 

Diagonalization have been performed for physically rel- 
evant values of A s _ of +220 K (= 150 cm -1 ) found for 
the Ti 3+ ion [2l| • The cubic CEF parameter B 4 is taken 
as +240 K (results are not sensitive to its exact value 
provided B4 > +50 K). The positive sign of B4 comes 
from ab initio point charge calculations of octupolar in- 
teractions of the Ti 3+ ion with the octahedral oxygen 
(negative charges) surroundings. Such value of B4 yields 
the T 2g -E g splitting of 2.5 eV. A splitting of 2.15-2.5 
eV has been observed for Ti 3+ ions in AI2O3, where the 
similar oxygen octahedron exists. In LaCoOa we derived 
value of B 4 of 280-320 K 

In Fig. 2 we show detailed calculations of the influ- 
ence of the spin-orbit coupling on the three lowest states 
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FIG. 1: The calculated fine electronic structure of the 3d 1 
electronic system (Ti 3+ , V 4+ ions) in the paramagnetic state 
under the action of the crystal field and spin-orbit interac- 
tions: a) the 10-fold degenerated 2 D term realized in the ab- 
sence of the CEF and the s-o interactions; b) the splitting of 
the 2 D term by the octahedral CEF surrounding E>4=+240 
K (A s _ =0) yielding the 2 T2 9 cubic subterm as the ground 
state; c) the splitting of the lowest 2 T2 S cubic subterm by the 
combined octahedral CEF and spin-orbit interactions (E>4= 
+240 K and A s _ = +220 K); the degeneracy and the asso- 
ciated magnetic moments are shown; d) the splitting due to 
the elongated tetragonal off-octahedral distortion of B|]= +10 
K (c/a>l); e) the splitting due to the compressing tetrago- 
nal distortion of B2 = -10 K (c/a<l, apical oxygens become 
closer). Figs c, d and e are not to the left hand energy scale 
- the splitting of the three lowest states on Figs c, d and e 
amounts to 333 K. 368 K and 372 K, respect ively. 



3d 1 




-""""^ B a = +240 K 






= +10 K 


■\ 




z = [001] 




3d' 




■^"^ B a = +240 K 
B° = -10 K 
z = [001] 


X 


~~~~2 — 
~~~~~2 





t- 



2y 



2 T2 g ) of the 3d 1 configuration for the elongated 



FIG. 2: The calculated spin-orbit coupling A dependence of 
the three lowest states t2 g (— 2 T2 g ) of the 3d 1 configuration 
for the elongated (c/a>l) and stretched (c/a<l) tetragonal 
off-octahedral distortion. The right hand states correspond 
to the three lowest states shown in Fig. Id and le. 

(c/a>l) and stretched (c/a<l) tetragonal off-octahedral 
distortion, realized for B° > and B° < 0, respectively. 
Analyzing the effect of the sign of the tetragonal distor- 
tion lead us to a conclusion that the magnetic moment 
lies along the tetragonal axis for the z-axis stretching. 
For the elongation case the moment lies perpendicularly 
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FIG. 3: The calculated temperature dependence of some 
properties of YTiOs. a) the temperature dependence of the 
three lowest states {ti g states) of the Ti 3+ -ion in YTiOs in 
the magnetically-ordered state below T c of 30.6 K; in the 
paramagnetic state the electronic structure is temperature 
independent unless we do consider a changing of the CEF 
parameters, for instance, due to the thermal lattice expan- 
sion. The used parameters: B 4 = +240 K, B" = -50 K, B|= 
-10 K, A s _ = +220 K and n d _ d = 69 T/fi B . The splitting 
of the Kramers doublets should be noticed in the ferromag- 
netic state. Excited states are at 377, 645, 28843 and 29447 
K. (b) the temperature dependence of the Ti 3+ -ion magnetic 
moment in YTiC>3. At K the total moment mn of 0.84 \ib 
is built up from the orbital m and spin m s moment of -0.06 
and 0.90 /j,b, respectively, c) The calculated temperature de- 
pendence of the 3d contribution Cd(T) to the heat capacity of 
YT1O3. The A- type peak marks T c . 



to the tetragonal axis. Taking into account that in YTiO,3 
the ordered moment lies along the c direction (in the 
Pbnm structure) and our long-lasting studies of CEF ef- 
fects [3 E3 113 HI we came to values for Bf] = -50 K 
and B 2 , = -10 K which reproduce the magnetic moment 
value, its direction and (some) thermodynamics. The 
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FIG. 4: Calculated temperature dependence of the paramag- 
netic susceptibility x (T) for the 3d 1 system Ti 3+ in YTi0 3 
for three different crystallographic directions (a, b, and c) cal- 
culated for B 4 = +240 K, B^= -50 K, B| = -10 K and A s _ 
= +220 K; these curves are denoted with CEF. The lowest 
solid line is the \c (T) dependence calculated with taking into 
account the ferromagnetic interactions with n d -d= 69 T//xb 
- this curve should be compared with experimental data. Ex- 
perimental data, taking after Ref. |4|, are shown by the low- 
est dashed line. The shadow area indicates the ferromagnetic 
state. The straight line denoted with S =1/2 shows the Curie 
law expected for the free S =1/2 spin. 

calculated ground-state eigenfunction (the z component 
of L and S are shown) 

iP GS ± = 0.690|±2,+i) - 0.678|+2,+i) 

- 0.253|±l,±i) - 0.020| + l,±i) (2) 

where the sign ± refers to two conjugate Kramers 
states. This state has S z =+0.44 and L z =± 0.10. The 
resultant moment m z =±0.78 fiB cancels each other in 
the paramagnetic state as is denoted in Fig. 3a. 

Making use of the \xy, +) function ex- 
tended for the spin component, as \xy,^f) — 
y/l/2 (|2,+^) — |— 2,+i)) (also functions \xz), 
\yz), \x 2 — y 2 ) and |z 2 ) ) one can write the ground state 
tpGS± function approximately as: 

iP GS ± = 0.967|xy,+) + 0.0085|x 2 -y 2 ,+) 

- 0.186|xz,±) - 0.172|?/z,±) + ....(3) 

In the magnetic state the ground state ipGS± Kramers 
doublet function becomes polarized as a molecular field 
is self-consistently settled down and the function 

"4>gs+ = 0.695|+2,-i) - 0.686|-2,-i) 

- 0.215| + l,+i) - 0.016|-l,+i) (4) 

is obtained as the ground state. The higher conjugate 
state is calculated to be described by 



-0G5_ = 0.665 1 
+ 0.3031—1, - 



1-2, +|) - 0.682| 
|) + 0.026|+1, 



-2,+i) 
~\) (5). 
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Below T c there opens, as is seen in Fig. 3a, a spin- 
like gap that amounts at T = K to 59.7 K (= 41.4 
cm -1 =5.15 meV). The spin- like gap is associated with 
the splitting of the Kramers doublet ground state in the 
ferromagnetic state. The magnetic ground state "0GS+ 
has S z = -0.45 and L z = +0.06 and the resultant mo- 
ment amounts to 0.84 \ib- The appearance of the mag- 
netic state is calculated self-consistently. It appears in 
the instability temperature (T c ) in the temperature de- 
pendence of the CEF paramagnetic susceptibility when 

Xcf( t c) = n d-d (6) 

as is illustrated in Fig. 4 for different crystallographic 
directions. With decreasing temperature this equality is 
reached the first for the c direction pointing the preferred 
magnetic arrangement axis. The magnetic state is calcu- 
lated self-consistently by adding to the Hamiltonian Eq. 
(1) the inter-site magnetic (spin-dependent) interactions 
instead of the last Zeeman term |29l | 

H d -d = n d -d {-m,d -m d + \ {m 2 d )) (7) 

where n d -d is the molecular-field coefficient. 

Having eigenvalues in a function of temperature we 
calculate the free energy F(T). From F(T) using well- 
known statistical formulae we calculate all thermody- 
namics like temperature dependence of the magnetic mo- 
ment, of the additional Cd heat capacity, of paramagnetic 
susceptibility, of the 3d-shell quadrupolar moment and 
many others. The present calculations are similar to 
those performed for FeBr2 30] and CoO 31]. 

CONCLUSIONS 

We have calculated consistently a value of the mag- 
netic moment of 0.84 fiB and its direction (along the c 
axis in the Pbnm structure) in YTiC>3 as well as temper- 
ature dependence of the paramagnetic susceptibility in 
very good agreement with experimental observations. 

This remarkable reproduction of so many physi- 
cal properties has been obtained within the localized- 
electron approach taking into account lattice off- 
octahedral distortions and the intra-atomic spin-orbit 
coupling [25j] . Although the spin-orbit coupling is weak, 
it amounts to only 1.2 % of the octahedral CEF splitting, 
it has enormous influence on the low-energy electronic 
structure and low-temperature magnetic and electronic 
properties. The spin-orbit coupling binds the orbital mo- 
ment to the spin moment being the reason for the spin- 
lattice and spin-phonon coupling. The suggested energy 
positions of the excited CEF states have to be verified by 
experiment. The good reproduction of the paramagnetic 
susceptibility indicates on the proper overall splitting of 
the tig states - one should note that our splitting is about 
5 times smaller than that obtained in Refs 0, 0, El • 
We are aware that our calculations have to be extended 
to take into account many other effects - first of them 
seems to be geometrical effects associated with the non- 
collinearity of the local symmetry axes, but these studies 



prove that magnetic properties of YTiC>3 are predomi- 
nantly determined by the atomic-scale lattice distortions, 
crystal-field and the spin-orbit coupling of the Ti 3+ ions, 
whereas charge fluctuations are negligible. An interplay 
of the spin-orbit coupling, lattice distortions and mag- 
netic order is very subtle involving rather small energies, 
smaller than 5 meV making theoretical studies quite dif- 
ficult. 

We point out that all discussed by us parameters are 
physical measurable parameters. The B° and Bj param- 
eters are related to the observed orthorhombic local sur- 
roundings. A negative value of the B!] parameter results 
from the stretching of the apical bonds with respect to 
the average bond length within the a ~ b plane. The B| 
parameter is related to the difference in the a— b plane. A 
success of our ionic approach, called due to extension to 
the magnetic state a quantum atomistic solid-state the- 
ory (QUASST) H3, if applied to YTi0 3 , is related to our 
long-lasting systematic studies, despite of discrimination 
and inquisition, of the spin-orbit coupling and crystal- 
field interactions in 3d compounds and in study of the 
region, where the spin-orbit coupling and off-octahedral 
lattice distortions are of the comparable strength. 
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